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Abstract. We show that the neutron star minimum pe- 
riod of uniform rotation is determined, for causal equa- 
tions of state, by the maximum value of the relativistic 
(compactness) parameter 2GM/ Rc 2 , allowed by causality 
for static neutron stars, 0.7081, and by the largest mea- 
sured mass of a neutron star. The relation between these 
three quantities, resulting from the extrapolation of the 
empirical formulae of Lasota et al. (1996), yields the min- 
imum period of 0.288 ms, only 2% higher than an absolute 
lower bound obtained in extensive exact numerical calcu- 
lations of Koranda et al. (1997). 

Key words: dense matter - stars: neutron - stars: pulsars 



1. Introduction 

The lower limit on the period of a uniformly rotating neu- 
tron star is sensitive to the equation of state (EOS) of 
dense matter above the nuclear density. Therefore, an un- 
certainty in the high density EOS implies a large uncer- 
tainty in the minimum period of uniform rotation, P m i n 
(see, e.g. Friedman & Ipser 1987; Friedman, Parker & Ipser 
1989; Salgado et al. 1994a,b; Cook et al. 1994). Hence, it is, 
therefore, of interest to find a lower limit on P m i n , that is 
independent of the EOS. This limit results from the con- 
dition of causality, combined with the requirement that 
EOS yields neutron stars with masses compatible with 
observed ones [currently the highest accurately measured 
neutron star mass is M"^ x = 1.442 M Q (Taylor & Weis- 
berg 1989)]. It will be hereafter referred to as P^in- 

The first calculation of -P^in was done by Glenden- 
ning (1992), who found the value of 0.33 ms. Glendenning 
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(1992), however, used a rather imprecise empirical for- 
mula, to calculate a lowest P m ; n by using the parameters 
(mass and radius) of the maximum mass configurations 
of a family of non-rotating neutron star models. His re- 
sult, therefore, should be considered only as an estimate of 
P^ n . Recently, Koranda et al. (1997) extracted the value 
of P^in from extensive exact calculations of uniformly ro- 
tating neutron star models. They have shown, that the 
method of Glendenning (1992) overestimated the value of 
P^ by 6%. The result of Koranda et al. (1997) calcula- 
tions can be summarized in a formula 

71 fmax 

P^k = 0.196 ms , 



(1) 



which combined with measured mass of PSR B1913+16 
yields today's lower bound for P^ = 0.282 ms. This ab- 
solute bound on the minimum period was obtained for the 
"causality limit (CL) EOS" p = (p — po)c 2 , which yields 
neutron star models of the surface density po and is max- 
imally stiff (dp/dp — c 2 ) everywhere within the star; it 
does not depend on the value of po- In the present letter 
we show that Eq. (1) can be reproduced using an empirical 
formula for P m i n derived for realistic causal EOS by La- 
sota et al. (1996), combined with an upper bound on the 
relativistic (compactness) parameter 2GM/Rc 2 for static 
neutron stars with causal EOS. 

2. Relation between x s and P m - m 

As shown by Lasota et al. (1996), numerical results of 
Salgado et al. (1994a, b) for the maximum frequency of 
uniform stable rotation can be reproduced (within better 
than 2%), for a broad set of realistic causal EOS of dense 
matter, by an empirical formula 
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where M s is the maximum mass of a spherical (nonrotat- 
ing) neutron star and R s is the corresponding radius, and 
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C{x s ) is a universal (i.e. independent of the EOS) func- 
tion of the compactness parameter x s = 2GM S /R s c 2 for 
the static maximum mass configuration, 



C(x s ) = 0.468 + 0.378x s . 

Combining Eq. (2) and Eq. (3) we get 



(P 



8.754 x 10~ 2 M, 



C(x s )x s 
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At given maximum mass of a spherical configuration, the 
maximum rotation frequency (minimum rotation period) 
is obtained for the maximum value of x s . At fixed x s , 
the value of P m i n is proportional to M s . Neutron stars 
for which masses have been measured, rotate so slowly 
that their structure can be very well approximated by that 
of a spherical star. Observations impose thus a condition 
M s > MX X - 

3. Lower bound on P m i U 

Our empirical relation, Eq. (4), indicates, that to minimize 
Pmin for given M"^ x we have to look for an EOS which 
yields maximum x s at M s = M^ x . It is well known, that 
if one relaxes the condition of causality, the absolute upper 
bound on x s for stable neutron star models is reached for 
an incompressible fluid (i.e., p = const.) EOS; the value 
of x a is then independent of M s and equal 8/9 (see, e.g., 
Shapiro & Teukolsky 1983). It is therefore rather natural 
to expect that in order to maximize x s under the condition 
of causality, one has to maximize sound velocity through- 
out the star. Together with condition of density continu- 
ity in the stellar interior this points out at the CL EOS, 
p = (p—po)c 2 , as to that which yields "maximally compact 
neutron stars"; introducing density discontinuities does 
not increase the value of Gondek & Zdunik (1995). 

[The conjecture that the CL EOS minimizes P m i n was al- 
ready proposed and then confirmed numerically in exten- 
sive exact calculations by Koranda et al. (1997)]. Note, 
that the value of x s for CL EOS does not depend on po 
(and therefore is M s -independent). It represents an abso- 
lute upper bound on x s for causal EOS, x s , max . Our nu- 
merical calculation gives x s (CL EOS) = a; s ,max = 0.7081. 
This corresponds to an absolute upper bound on the sur- 
face redshift of neutron star models with causal EOS, 
W = (1 - a^max)- 1/2 - 1 = 0.8509. 

Let us consider the effect of the presence of a crust 
(more generally, of an envelope of normal neutron star 
matter). For a given EOS of the normal envelope, the rel- 
evant (small) parameter is the ratio pb/pbC 2 , where pb and 
Pb are, respectively, pressure and mass density at the bot- 
tom of the crust (Lindblom 1984). The case of pb = cor- 
responds to stellar models with no normal crust. Numeri- 
cal calculations show, that adding a crust onto a CL EOS 
core implies an increase of R s , which is linear in pb/pbC 2 ; 



for a solid crust we have typically Pb/pbC 2 ~ 10~ 2 . The 
change (increase) in M s is negligibly small; it turns out to 
be quadratic in Pb/pbC 2 . This implies, that the decrease of 



and of the maximum surface redshift z s 



due to 



the presence of a crust, is proportional to Pb/pbC 2 . This is 
consistent with Table 1 of Lindblom (1984). However, the 
extrapolation of his results to pb = yields z max = 0.891, 
which is nearly 5% higher than our value of z max ! This 
might reflect a lack of precision of the variational method 
used by Lindblom (1984), which led to an overestimate 
of the value of z max . It should be stressed that while a 
precise determination of M max = M s for static neutron 
star models is rather easy, determination of the precise 
value value of the radius of the maximum mass configu- 
ration, i? s , (with the same relative precision as M s ) and 
consequently of the value of x s (with, say, four significant 
digits), is much more difficult and requires a rather high 
precision of numerical integration of the TOV equations. 

In what follows, we restrict ourselves to the case of the 
absolute upper bound on x s , obtained for neutron star 
models with no crust. Inserting the value of x Sjmax into 
Eq. (4) we get 
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Current lower bound on P, resulting from the above equa- 
tion, is thus 0.288 ms, which is only 2% higher than the re- 
sult of extensive exact numerical calculations of Koranda 
et al. (1997). 

The formula (^|) deserves an additional comment. In 
numerical calculations, of a family of stable uniformly ro- 
tating stellar models, for a given EOS of dense matter, 
one has to distinguish between the rotating configura- 
tion of maximum mass, which corresponds to the rotation 
frequency fW max (EOS), and the maximally rotating one, 
which rotates at f2 max (EOS) (Cook et al. 1994, Stergioulas 
& Friedman 1995). Notice, that determination of a maxi- 
mum mass rotating configuration (and therefore of f2A/ max ) 
is a much simpler task than the calculation of exact value 
of f2 max , which is time consuming and very demanding as 
far as the precision of numerical calculations is concerned. 
Usually, both configurations are very close to each other, 
and O max is typically only 1-2% higher than fW max ; such 
a small difference is within the typical precision of the em- 
pirical formulae for !l m „. Actually, the formula for C(x s ), 
Eq. (3), was fitted to the values of f^/ max (EOS) calcu- 
lated in (Salgado et al. 1994a,b). Therefore, Eq.(||) should 
in principle be used to evaluate the causal lower bound 
to P m in,A/ max ; it actually reproduces, within 0.2%, the ex- 
act formula for this quantity, obtained by Koranda et al. 
(1997) [see their Eq. (8)]. 

It should be stressed that Eq. (5) results from an ex- 
trapolation of the empirical formula of Lasota et al. (1996). 
General experience shows that - in contrast to interpola- 
tion - extrapolation is a risky procedure. The fact that 
in our case extrapolation of an empirical formula yields - 
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within 2% - the value of P m in of Koranda et al. (1997) (and 
reproduces their value of P m in,M max ), proves the usefulness 
of compact "empirical expressions" which might summa- 
rize, in a quantitative way, a relevant content of extensive 
numerical calculations of uniformly rotating neutron star 
models. 
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